Effect of thermal noise on vesicles and capsules in shear flow 
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We add thermal noise consistently to reduced models of undeformable vesicles and capsules in 
shear flow and derive analytically the corresponding stochastic equations of motion. We calculate 
the steady-state probability distribution function and construct the corresponding phase diagrams 
for the different dynamical regimes. For fluid vesicles, we predict that at small shear rates thermal 
fluctuations induce a tumbling motion for any viscosity contrast. For elastic capsules, due to thermal 
mixing, an intermittent regime appears in regions where deterministic models predict only pure tank 
treading or tumbling. 
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Introduction.- The dynamics of fluid vesicles and elas- 
tic capsules in linear shear flow has been extensively stud- 
ied in the past decades (for reviews, see Refs. [IHI])- For 
fluid vesicles, the transition between the tank-treading 
regime (TT), in which the membrane rotates around the 
interior fluid at constant orientation, and the tumbling 
regime (TB), where the capsule rotates as a rigid body, 
is described with relatively good accuracy by the Keller- 
Skalak (KS) model [5], which assumes ellipsoidal vesi- 
cles of fixed shape. This assumption holds for moderate 
shear rates, while other dynamical regimes are observed 
if the vesicle deforms due to higher shear stresses [6]- 
[15] • The Skotheim-Secomb (SS) model [16] extends the 
KS model to capsules showing elastic behavior, such as 
the shape memory of red blood cells [17 . It captures 
many features observed in experiments [18] and simula- 
tions [T9H22] . For capsules, the TT motion exhibits a pe- 
riodic oscillation of the orientation angle called swinging 
(SW). The TB regime also shows a periodic oscillation 
of the tank-treading angle. In addition, an intermittent 
regime (INT) appears for which both swinging and tum- 
bling happen alternately. A thorough analysis of this 
intermittent regime is found in Ref. [23] , which suggests 
that intermittency should also exist for deformable cap- 
sules. However, analytical work [4j |24] and simulations 
[T9H2T| [25] tend to show that the intermittent regime 
may be only transient if capsule deformation is consis- 
tently taken into account. 

Since these objects are micrometer sized, it can be rea- 
sonably assumed that thermal fluctuations should have 
an influence on the dynamical transitions. Indeed, ex- 
periments on vesicles [26j [27] show that they lead to 
sensible discrepancies from deterministic theoretical de- 
scriptions. However, most theoretical work that includes 
thermal noise has been done either for quasi-spherical 
vesicles in two dimensions (2D) [28] [29] and three di- 
mensions (3D) [30] or through pure phenomenological 
models [9j [31j [32]. Thermodynamic consistency, how- 
ever, imposes strong conditions on the stochastic equa- 
tions of motion which we derive here in a minimal ex- 
tension of the SS model. We then determine how such 
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FIG. 1. Schematic view of a capsule in shear flow. 



thermal fluctuations affect the phase diagram for vesicles 
and capsules. In particular, we show that under realistic 
conditions thermal effects should be observable. 

Reduced model- We describe the dynamics of an un- 
deformable, ellipsoidal vesicle or capsule in linear shear 
flow as represented in Fig. [T] The two degrees of free- 
dom are the inclination angle of the particle relative 
to the middle plane of the shear flow and the phase an- 
gle (j) describing the tank-treading motion of the outer 
membrane. Following Ref. [16] . the balances of angular 
momentum and dissipated energy lead to the equations 
of motion 

d t = -l(l-^co S (20))-^+Cs(t), (1) 
2 \ z J z 

W = ~4tJ (7- cos(20) + -% sin(20)) + 

(2) 

which give the temporal evolution of and (j) as functions 
of the shear rate 7, the viscosity ratio A = r]i/r] between 
inner and outer viscosity, the main axes ai, <22, a 3 and 
the volume V of the capsule, the elastic energy Eq cor- 
responding to the shape- memory effect, and geometrical 
quantities [33] . We quantify the influence of thermal fluc- 
tuations by adding the stochastic noises (o(t) and C<f>(t) 
similarly to Refs. [28] [31]. They obey 

<Ci(*)> = 0, <C*(W)> = 2AM(* - 0, (3) 
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where {i,j} G 0} and Dq and are the respective 
diffusion constants. The rotational diffusion constant is 
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which follows from the response of a rigid ellipsoid to a 
torque about the z axis in a liquid of viscosity r\ [34] , k B 
being the Boltzmann constant, T the temperature, and 
diid2 geometrical quantities [35]. To determine D^, we 
use the fact that the shape-memory effect was introduced 
by adding the potential Eq sin 2 <\> to the energy balance 
of the capsule. The fluctuation-dissipation relation in 
equilibrium (7 = 0) applied to Eq. ([2| then implies 
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The diffusion coefficients Q and ^ depend on the in- 
verse of the capsule volume V, i. e., the effects of thermal 
fluctuations become more important when the capsules 
are smaller. A small outer viscosity r] also amplifies 
thermal effects. 

Equations ([I]) and Q form a set of Langevin equations 
which can be transformed into a single Fokker-Planck 
equation for the probability distribution p(9,(j),t) [36] . 
This Fokker-Planck equation reads 



d t p(0, <j>, t) = -d e je - <9 j , 
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with the currents 
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where we defined B = zi/zq + 2/\'zqZiZ2 and 
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In the following, we will solve Eq. (|6| in the steady state, 
i. e., for d t p(0, 4>, t) = 0. 

In the absence of shear flow (7 = 0), the equilibrium 
probability distribution becomes 
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where Iq(x) = 1/ir J Q exp(xcos^) dO is a modified Bessel 
function of the first kind. This distribution does not de- 
pend on since no orientation is preferred. For fluid 



inelastic vesicles we have Eq = 0, leading to a homoge- 
neous distribution for both angles. For elastic capsules, 
the distribution shows two peaks at the angles = and 
(f) = 7r due to the shape-memory energy Eq sin 2 0. 

In order to investigate the nonequilibrium steady states 
(7 > 0) and construct phase diagrams, we introduce the 
tumbling ratio [37] 
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which plays the role of an order parameter: uj tu = 
corresponds to a pure tank-treading motion and uj tu — 1 
to pure tumbling, intermediate values corresponding to a 
dynamics in which both tumbling and tank treading take 
place. In the following, all numerical calculations are 
performed by setting a\ = as = 4 um and CL2 = 1.5 um, 
which corresponds approximately to the dimensions of 
a red blood cell [18 . The outer viscosity is taken to 
be approximately that of water, r] = ImPas, and the 
temperature T = 293 K. 

Fluid vesicles - Fluid vesicles do not exhibit shape 
memory, i. e., Eq = 0, implying that the Fokker-Planck 
equation ([6| does not contain any term depending explic- 
itly on (j). Therefore, we can write p s {0 1 (f)) = p s (0). In 
addition, by averaging Eqs. ([T]) and one obtains 
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such that the tank-treading frequency {d t <j)) can be de- 
duced from the tumbling frequency (dfO). 

Without thermal fluctuations, we recover the KS 
model for which two cases must be distinguished. If 
B < 1, i. e., above the critical viscosity contrast 
A c = 1 + 2zq jz\ Z2, the solution of Eq. ^ is 
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(13) 



which corresponds to a tumbling motion (TB) of the cap- 
sule around the z axis. The mean tumbling frequency 
calculated from ([T]) is equal to (d t 0) = — 7VI — B 2 /2. 
Note that even though the capsule tumbles, the mean 



tank-treading velocity ( 12 ) becomes zero only in the limit 
A — >• 00, a fact already stressed in Ref. [5 . The TB 
regime for vesicles therefore consists in a mixing between 
tumbling and tank treading. If B > 1, i. e., A < A c , the 
inclination angle takes the stationary value 
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while the membrane exhibits a tank-treading motion 
(TT) around the interior of the capsule. Here we have 
(d t 0) = and (dt4>) = — A f/\ f Z\Z2B, which corresponds 
indeed to pure tank treading without tumbling. This 
transition between tumbling and tank treading happens 
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FIG. 2. Phase diagram for the KS model with thermal noise 
as a function of 7 and A. Inset: cj tu along the white dotted 
line (7 = Is -1 ), and the dashed line corresponds to the deter- 
ministic KS model. An enlargement of the bottom left corner 
is also shown. 



at the critical viscosity ratio A c which depends only on 
the geometry of the vesicle but not on the shear rate 
7. The order parameter (11) is u tu = for TT, and 
< cu tu < 1 for TB. 



If we include thermal noise, the Fokker-Planck equa- 
tion (pj can still be easily solved as 



where 
J\f = 27r / exp 
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is the normalization constant and Iq the Bessel function 



defined after Eq. (10). In this case, the stationary proba- 
bility distribution ( 15 ) depends on the shear rate 7, which 
was not the case without noise. The mean tumbling fre- 
quency 
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is never exactly zero, such that co tu ^ always. How- 
ever, as can be seen in the inset of Fig. [2] showing the 
tumbling ratio as a function of A for 7 = 1 s -1 , it decays 
exponentially in the black region. We can thus still effec- 
tively identify a TT regime, defined here by uj t u < 0.1, in 
which tumbling becomes rare. The contour plot of Fig. 
[2| which represents u tu as a function of A and 7, can 
therefore be interpreted as a phase diagram, where TT 
occurs in the black region and TB elsewhere. 

The first effect of thermal noise is to smooth out the 
TT-TB transition, which happens for smaller viscosity 
contrasts than in the deterministic case, as can be seen in 
the inset of Fig. [2] Moreover, in this case, this transition 
depends not only on A but also on 7. In the limit of high 
shear rates, fluctuations are suppressed and we recover 
the critical viscosity contrast of the KS model. On the 
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FIG. 3. (Color online) (a) Probability distribution of the ori- 
entation angle at A = 2 and different shear rates, (b) Mean 
inclination angle (0) as a function of the viscosity contrast. 
The solid black lines represent the KS model. 



other hand, the critical viscosity contrast goes to for 
low but non-zero shear rates. In fact, as the enlargement 
in Fig. [2] shows, a TB motion is always present at shear 
rates lower than approximately 0.1s -1 , which is a new 
feature with respect to the deterministic phase diagram 
derived in Ref. [TO] . 

In order to illustrate this shear-induced phase transi- 
tion, we plot the distribution function p s (0) at A = 2 and 
for 7 = 0.1, 1, and 3 s -1 in Fig. [3|a). In the determin- 
istic model, the vesicle tank treads at fixed orientation 
#tt (vertical line). With thermal noise, the dynamics is 
different for these three shear rates indicated by crosses 
in Fig. [2] For 7 = 0.1s -1 , the vesicle's orientation can- 
not remain fixed because most fluctuations induce ro- 
tations and the regime is effectively TB (u tu = 0.68). 
The distribution is then relatively flat (dashed-dotted 
line). For 7 = Is -1 , the vesicle is tank treading and 
from time to time a full rotation is induced by fluctu- 
ations (cjtu = 0.13). The probability distribution is al- 
most symmetrically centered around #tt (dashed line). 
For 7 = 3s -1 , the regime is effectively TT since tum- 
bling practically never happens (uj tu = 0.0052) and the 
orientation angle fluctuates around #tt due to thermal 
noise (dotted line). The mean inclination angle (6) as a 
function of A is also influenced by thermal noise - a fact 
already seen in Ref. [31] - and an exponential tail appears 
instead of sharp transition, as illustrated by Fig. |3|b). 
Such a tail is similar to the one observed in a simulation 
including thermal noise [29]. A similar slow decay of the 
mean inclination angle has also been recently observed 
in an experiment [6j [27] , but it seems that it can also be 
reproduced in deterministic simulations [15]. 

Elastic capsules - For Eo > 0, the solution of (|6| can- 
not be separated and one has to solve the equation nu- 
merically. The phase diagrams of Fig. [4] represent the 
tumbling ratio uj t u as a function of the viscosity contrast 
A and the shape-memory energy Eq at a shear rate of 
7 = Is -1 . As for vesicles, u tu is never exactly or 1 
with thermal noise, but the values are extremely close 
in the black and white regions such that we can still ef- 
fectively identify the SW and TB regimes, respectively. 
The region for which intermittent dynamics takes place 



FIG. 4. (Color online) Phase diagrams of the dynamics of 
capsules as a function of A and Eo (in units of k^T) for 7 = 
Is -1 for the deterministic (left) and stochastic (right) case. 




FIG. 5. (Color online) (a) cj t u as a function of Eo. (b)-(d) 
Stationary probability distribution p s (6, <j)) for three different 
values of Eq. The other parameters are A = and 7 = 1 s _1 . 



is much wider in the presence of thermal noise. In par- 
ticular, there is now a large intermittent region for A < 1 
which was not there in the deterministic case. 

This can be seen in Fig. [5j which shows the phase 
transition for A = 0. In the SS model, this transition is 
discontinuous and happens between Eo = 96 (in units 
of k&T), where the capsule is perfectly tank treading 
(cjtu = 1)? an d Eo = 97, at which perfect tumbling 
(cjtu = 0) is observed [see Fig. (5ja)]. The thermal 
fluctuations turn the transition into a continuous one, 
which implies a mixing of both types of motion in a re- 
gion of approximately 40/cbT around the critical energy. 
This is illustrated by the stationary probability distribu- 
tion at different values of the shape-memory energy. For 
Eq = 60 [Fig. [5jb)], the capsule is in the TT regime 
and the probability distribution follows exactly the de- 
terministic trajectory (dashed line), but the amplitude 
of the oscillations in 6 is approximately two times larger. 
Similarly, for Eq = 120 [Fig. [5jc)], the capsule is tank 
treading and the most probable path is the determinis- 
tic one (dotted line) with larger oscillations of (j). At the 
transition [Eq = 96.5, Fig. [5]i], however, we observe a 



FIG. 6. (Color online) Intermittent behavior with 7 = Is -1 , 
A = 4 and Eo — 40/cbT. Left: Deterministic trajectory of 
length 500 s. Right: Corresponding probability distribution 
in the presence of thermal noise. 



thermal mixing of the TT and TB regimes. The capsule 
exhibits both motions with a slight offset with respect 
to the deterministic paths at Eq = 96 (dashed line) and 
E = 97 (dotted line). 

In the presence of thermal noise, the intermittent 
regime predicted by the SS model thus still exists, but 
is superimposed with a thermal mixing of the dynamical 
regimes, as shown in Fig. [6j The left panel shows one 
deterministic trajectory of the capsule in the — (j) phase 
space. The trajectory is quasi-periodic and consists of 
incomplete tumbling rotations (horizontal curves) inter- 
rupted by incomplete tank-treading motions since there 
must not be any intersection. The right panel shows the 
probability distribution in the presence of thermal noise. 
The tank-treading path is the most probable one, even 
though tumbling also happens. The white areas on the 
top right and top left are almost excluded, which was 
not the case for the deterministic trajectory. Complete 
TB rotations as well as complete TT rotations may oc- 
cur since an intersecting trajectory is not forbidden for a 
stochastic motion. 

Concluding Perspectives.- We have shown that under 
realistic conditions, the effects of thermal fluctuations on 
the dynamics of fluid vesicles and elastic capsules should 
be observable for objects having approximately the same 
dimensions as red blood cells in water. In most extant ex- 
periments on elastic capsules, either the capsules are too 
big [38] or the outside medium has a too high viscosity 
[T8] [39] for showing thermal effects. Two important pre- 
dictions concern the tumbling motion of fluid vesicles at 
low shear rates for all viscosity contrasts, and the thermal 
mixing of swinging and tumbling for elastic capsules at 
low viscosity contrast. In the future, the generalization of 
our model to deformable objects should be investigated, 
for which we expect that thermal noise would play an 
even more important role, especially for fluid vesicles as 
suggested by the experiments reported in Ref. [27] . 
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